Coexistence of vector chiral order and Tomonaga-Luttinger liquid in the frustrated 

three-leg spin tube in a magnetic field 
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The frustrated three- leg antiferromagnetic spin-^ tube with a weak interchain couphng in a mag- 
netic field is investigated by means of Abelian bosonization techniques. It is clearly shown that a 
vector chiral long-range order and a one-component Tomonaga-Luttinger liquid coexist in a wide 
magnetic-field region from a state with a small magnetization to a nearly saturated one. The chiral 
order is predicted to still survive in the intermediate plateau state. We further predict that (even) 
when the strength of one bond in the three rung couplings is decreased (increased), an Ising-type 
quantum phase transition takes place and the chirality vanishes (no singular phenomena occur and 
the chiral order is maintained). Even without magnetic fields, the chiral order would also be present 
if the spin tube possesses easy-plane anisotropy. 

PACS numbers; 75.10.Jm, 75.40.Cx, 75.50.Ee 



I. INTRODUCTION 

Quantum spin tubes— namely, systems of coupled spin 
chains with a periodic boundary condition (PBC) for the 
interchain (rung) direction— have been studied since late 
1990sJ-^i^iii^i^:^iSaaai in particular, odd-leg spin tubes 
have attracted much interest, because the antiferromag- 
netic (AF) rung coupling yields frustration. It has been 
recognized now that coupled spin chains including spin 
tubes offer fairly richer physics compared with single spin 
chainsJ^ Recently, a few new spin-tube compounds such 
as a three-leg tube^ [(CuCl2tachII)3Cl]Cl2 and a nine- 
leg one^ Na2V307 have been fabricated. The presence of 
these materials has further promoted the study of spin 
tubes. 

For nonfrustrated AF spin-S" ladders— namely, cou- 
pled AF chains with an open boundary condition for 
the rung— the following universal feature (generaliza- 
tion of Haldane's prediction) has been elucidated i^^'^"^i^^ 
if both 25* and the number of chains, N, are odd, 
the low-energy physics is described by a gapless, one- 
component Tomonaga-Luttinger liquid (TLL) and no 
symmetry breakings take place, whereas for all the other 
cases (25* xN — even) , the system has a gapped spectrum 
with conserving all symmetries. This even-odd property, 
however, does not always hold in frustrated coupled spin 
chains. In the three-leg AF spin-^ tube where the rung 
triangle induces geometrical frustration, it is knowi*^ that 
the one-site translational symmetry along the chain di- 
rection is spontaneously broken and a finite gap between 
a ground state and the lowest excitation exists. 

In this paper, we point out another striking difference 
between the three-leg nonfrustrated ladder and the frus- 
trated tube: namely, utilizing bosonization techniques, 
we show in a clear way that a vector chiral long-range 
order resides in the magnetic-field-induced TLL phase 
of the three-leg AF spin-i tube at least in the weak- 
rung-coupling regime. In the chiral phase, the parity 
symmetry for the rung direction is spontaneously bro- 



ken. On the contrary, the chiral order is absent in 
the corresponding ladder. For one-dimensional gapped 
U (l)-symmetric coupled spin chains, (i) the field- induced 
magnon-condensed gapless state is usually described by a 
standard TLL and conserves all symmetries and (ii) many 
physicists have focused only on whether nontrivial mag- 
netization plateau states exist or not. Therefore, the pos- 
sibility of chirality in spin tubes has not been considered 
well so far.™ We also show that the chiral order continues 
in the intermediate plateau state of the spin tube. These 
predictions means that the three-leg spin tube possesses 
both gapless and gapped chiral-ordered states. We fur- 
ther discuss a rung-coupling modification in the three-leg 
spin tube: the strength of one bond of three rung cou- 
plings is changed. If the strength of the changed bond is 
extremely increased (decreased), the system approaches 
a two-leg ladder plus single chain (a three-leg ladder). 
Both limits possess no geometrical frustration. Interest- 
ingly, it is found that in the case of increasing the bond 
strength, the chiral order always remains unbroken, but 
an Ising-type quantum phase transition takes place and 
the chirality disappears at the Ising critical point in the 
other case. 

The organization of the paper is as follows. In Sec. [TTl 
we define a spin tube model and briefly explain the 
bosonization method. The next two sections are the main 
part of the paper. We show the mechanism of the vector 
chiral order in Sec. IIIII Section IIVI is devoted to inves- 
tigating effects of the rung deformation. Finally, conclu- 
sions and brief discussions about our results are presented 
in Sec. El 



II. 



MODEL AND METHOD 



In this paper, we mainly consider the three-leg AF 
spin-i tube with a homogeneous rung coupling. The 



2 



Hamiltonian is defined as 



Ti. = l^'^'Snj' ■ Sn,j + 1 + JA_Sn,j ' <5Vi+lj " HS^ j ,(1) 



where Snj is spin-i operator on site j of the nth chain 
[n = 1-3), J > (J_L > 0) is the AF intrachain 
(rung) coupHng, and the PBC S4J = Sij is imposed. 
We begin with the independent three chains in a field 
H. From Abehan bosonization^i^ the nth chain in the 
low-energy limit is mapped to a Gaussian model 7i„ ~ 
J dx ^ [K{dx0n)^ + K~^{dx4>n)'^] , where (/)„ is the scalar 
boson field and On is the dual to (f>n- The TLL param- 
eter K and the spin-wave velocity v depend on J and 
H; when H is changed fi-om to the saturation field 2 J, 
V {K) monotonically decreases (increases) from 7rJao/2 
(1/2) to (1) [ao- lattice constant]. The spin operators 
are also bosonized as 



M + aodx(t>n/\^ 



+ {-iyAi sin(V47r0„ -I- 27rAfj) -I- 



S, 



Bi sm{V4:TT(f)n + 27rA/j) 



(2a) 



(2b) 



where M{H) = {S^^), and Ai and Bi are nonuniver- 
sal constants of 0(l)fiiii^ This formula indicates that 
the period of (pn {On) is -x/tt (-s/iTr). Using 7i„ and for- 
mula ([2]), we can straightforwardly derive the following 
effective Hamiltonian of the model ([T]) in the weak-rung- 
coupling regime: 

f 3 

T^eff - dxJ2l [K{dxOn)^+K-\dxCj)n?] 

^ 2MJ_L Q^^^ JlOo y-^ dx4'ndx<j),i+l 
^ n n 

+ XI '^^^^^ COs(V7r(6l„ - ein+i)) 



+ cos(\/47r(0„ - 0,i+i)) 



(3) 



where (^4 = 0i, ^4 = 0i, and we have written only the 
important part among all the rung-coupling terms and 
neglected terms with oscillating factors e'*^^^-* or (— l)-*, 
which are irrelevant. In the bosonization picture, sym- 
metries of the spin tube H]) are represented as follows: 
a J7(l) rotation around the S"^ axis S^,j e^'^S^j, the 
one-site translation along the chain j Sn that 
along the rung — > S"_^_i j, the site-parity transfor- 
mation for the chain S" j — > S"^_j, and that for the 
rung j ^ 5*3 J, respectively, correspond to On 

On + i<Pn, On) ^ (</>„ + V^(M + 1/2), 6I„ + 0F), 

((/>„, ^n) {(l)n+l,On+l), {<pn{x) , On{x)) (-(/i„(-x) + 

y/Tr/2,0n{-x)) and (01,6*1) ^ {(1)3, Ss)- These sym- 
metries strongly restrict possible terms in Hcs- In 




FIG. 1: Relationship between {9n} and {O5}. C is a constant. 



other words, the symmetries reduce the number of cou- 
pling constants. As a result, for all vertex opera- 
tors without oscillating factors and conformal spin, only 
X;„cos[/V47r((/>« - (/>«+i)] and cos[/\/47r(6'„ - On+i)] 
are allowed to exist in Eq. ([3]). The most relevant terms 
with I = 1 indeed appear in Eq. ([3]). Utilizing the ef- 
fective theory we study the low-energy properties 
of the spin tube below. We stress here that almost 
all the discussions below are applicable to other three- 
leg spin tubes with symmetry-conserved perturbations 
{XXZ anisotropy, further-neighbor interactions, etc.). 



III. 



VECTOR CHIRAL ORDER 



Let us introduce new boson fields <I>o — '^n^^/V^, 

$1 = (01 - 03)/%/2, $2 = (01 + 03 - 202)/v^, 60 = 

E„^n/V3, 61 = {Oi - 03)/V2, and 62 = (^i + ^3 - 
26*2)/ V6- The relationship between old and new bosons is 
illustrated in Fig.[TJ As one will see later, these fields help 
us to detect a vector chiral order. Using these new fields, 
we can diagonalize the boson bilinear part in Eq. ([3]). 
Consequently, the effective Hamiltonian is rewritten as 

r ^ 

Heft ^ dx [Kg{dxeg)' + K-\dx<^gf] 

q=Q 



12 J± 



+Af-^y[$l,$2] + -- - , 

2ao 

where the potential I/[q;, /3] is defined as 
y = 2cos ^ 



(4) 



M|a)cos(^^/3) -Hcos(\/2^a) . (5) 



The new TLL parameters Kq and velocities Vq are eval- 
uated as 



Kq = KJq \ Vq = vfo, 
Kl^2 = KJq^ = Kg, Vi:2 = vfg = V, 



9' 



(6) 



where /o = (1 + 2f ^)i/2 and fg^{l- 

One finds that the Hamiltonian (|3|) does not contain 
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FIG. 2: (a) Potential V[Qi,e2]- (b) Physically relevant zone 
(enclosed by the green line) in the projected Q1-Q2 plane. 
Kn play the role of the primitive translation vectors. The 
rung-parity transformation causes Oi — 6i, and the rung 
translation by one site does K„ Kn+i- 



any vertex operator with $0 or Qq. This is strongly 
supported by symmetries of the tube. The $0 sector 
hence provides a one-component TLL. The derivative 
term dx^o can be absorbed into the Gaussian part and 
yields only a small correction of M. On the other hand, 
the remaining $1^2 sectors are subject to the effect of the 
potential V. Since the scaling dimensions of V^[<i>i,$2] 
and ^[61,92] are, respectively, evaluated as 2Kg and 
l/{2Kg), the latter F[0i,62] is always dominant and 
relevant [i.e., l/{2Kg) < 2] from a small M to the satu- 
ration. The phase fields Qi and 82 thus are pinned at 
the minimum of y [01,02], and $12 sectors have mas- 
sive spectra. Consequently, the low-energy physics of the 
tube (dJ is governed by a TLL in the $0 sector. The 
TLL just corresponds to the known field-induced critical 
phase.^ [For the case of M = 0, e''^*^^ = 1 and rele- 
vant vertex operators with $0 appear in Eq. They 
are expected to cause the $0 sector to obtain a gapped 
spectrum.] 

The scenario leading to the TLL here is already known 
well^2,i2, However, as one will see below, characteristics 
of the frustrated spin tube ([1]) are hidden in the form 
of the potential V. Figure ^a) and Eq. ^ show that 
six minimum points {Qi, Q2) — (±V27r/3, ±-y/27r/3) and 
(±2-\/27r/3, 0) lie around the origin (0,0). [Inversely, 
(0,0) is the unique minimum of the potential J±V for 
Jj_ < 0.] Not all, however, are physically meaningful, be- 
cause phase fields 0„ have a period. Projecting the phys- 
ically relevant cubic space 0„ € [— -^7?, (mod ^/4tt) 
onto the 81-62 plane (see Fig. [1]), we obtain a diamond 
zone as in Fig. [2jb). Therefore, it is enough to consider 
only two minimum points (8^,82) = (±\/2^/3, ^27r/3) 
that satisfy 8^ — —8^ ^ 0. Since the rung-parity trans- 
formation {<j)i,6i) ^ (03,^3) causes 81 —81, pinning 
81 to these two minima implies the spontaneous break- 
down of the rung-parity symmetry. As a candidate of 
the rung-parity order parameter, we can propose a vec- 
tor chirality k„j = Sn,j x S„+ij, which changes sign via 
the rung-parity operation. From the formula ([2]), the z 
component of a chirality is evaluated as 

{kI^) « {|sin (%/2^8i) ^ + • • • = T finite, (7) 



for 8*. Remarkably, the leading term of the chiral- 
ity does not contain the massless fields ($0)8o). Sim- 
ilarly, {n\(2) j) shown to be equivalent to (n^j)- 
We thus conclude that the vector chiral long-range or- 
der exists in the field-induced TLL phase. The chiral- 
ity correlation function exhibits an exponential decay: 

(«3,,«3,o) « Kj)^ + C^oe-I^l/VTM + • • ■ ■ Since any 
artificial approximation such as a mean-field decoupling 
is not used to obtain Eq. ([7]) , the prediction of the chiral 
order is highly reliable. Note that in the chiral phase, 
the spontaneous breakdown of the Z2 rung-parity sym- 
metry occurs, but the U{1) spin-rotational symmetry is 
preserved. 

For the spin-i tube ([1]), it is predicted that the inter- 
mediate plateau state appears at M = 1/6 in the region 
Jl > Jli^ 0.1J);i In the bosonization picture, the spin 
gap of the plateau state is interpreted as that of the $0 
sector, which is induced by a relevant higher-order umk- 
lapp term cos(4\/37r$o + 127rMj) = cos(4\/37r$o)- On 
the other hand, the $1^2 sectors are hardly influenced by 
whether the $0 sector obtains a gapped spectrum or not. 
Therefore, we predict that the chiral order still continues 
in the intermediate plateau phase at least if Jl(> J^) 
is sufiiciently small. This chiral plateau state is reminis- 
cent of the narrow chiral phase in the classical AF XY 
model on a triangular lattice^^ It is noteworthy that by 
controlling the strength of the magnetic field, one can re- 
alize both gapless and gapped chiral states in single spin 
tube. 



IV. RUNG DEFORMATION 

In this section, we discuss the rung deformation. Let 
us modify the coupling between the first and third chains 
as 

J±Si.rS3,j -> J^{l + 5)Si.j-S3j. (8) 

This explicitly violates the rung translational symme- 
try, but conserves the parity symmetry between the first 
and third chains. Therefore, is still valid as a 

rung-parity order parameter. The cases oi 6 — —I 
and 6 +00, respectively, correspond to a three- leg 
ladder and a system of a two-leg ladder plus single 
chain. Note here that if we set J±S = const <C J, the 
present weak-rung-coupling approach is available even 
for the case of a large |(5| 3> 1. A finite S brings new 
bosonic terms 5J±Mdx{<pi+<t'3)/\/7T, SJ±aodx4>idx(l)3/Tr, 
(5Jj_-B§ao cos('\/27r8i), etc. The first term can be ab- 
sorbed into the Gaussian part via the shift of which 
does not affect dual fields 8, and yields a small deviation 
from the uniform magnetization, {S2 j) ^ {^1(3) j)- I* 
not surprising because (as mentioned above) S breaks the 
rung translational symmetry. The second boson bilinear 
term is expected not to qualitatively influence the low- 
energy physics. Actually, introducing a new basis dif- 
ferent from (<I>g,8g), we can diagonalize all the bilinear 
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FIG. 3: Deformed potential Fa [61,82] 



terms. The main effect of the rung modification origi- 
nates from the third cosine term which varies the form 
of the potential T^[6i, O2] as follows: 



V[Qi,e2] ^ V5[Qi,Q2] = 



(9) 



2 cos 



81 ] cos 



(1 + 6) cos V27rei 



The deformed potential Vs is illustrated in Fig. [31 From 
Vs, one finds that when S increases [decreases and reaches 
del = ^0.5], the two minimum points still survive with 
keeping Oi 7^ [approach with each other and meet at 
(61,62) = (0, y/2Tr/3)]. The chirality is hence 

shown to be always finite when an arbitrary positive S is 
applied. In other words, once an infinitesimal exchange 
is added between a two-leg AF ladder and single chain, 
a vector chiral order immediately emerges. 6n the other 
hand, for the case of 5 < 0, the vector chiral order is 
predicted to vanish at a certain point 6 = Sc and the low- 
energy physics is described by the TLL of the $0 sector 
for S < 6c- The true value of 6c must be renormalized 
from the classical value 6c\ by the effect of the Gaussian 
part and must depend on J±/J and H/J. 

In order to further understand the physics near the 
phase transition ai d — 6c, it is necessary to go beyond 
the semiclassical analysis of Vs. If cos(\/37r/262) in Vs 
is allowed to be replaced with its expectation value, the 
effective Hamiltonian for the $1 sector is written as the 
following double sine-Gordon (dSG) mode&^: 

H[*,,ei] ~ Jdx [Kgid^Q.f + K;\d.,'Pif] 



Ci cos (W |©i) - (1 + ^) COS (v^ei) 



(10) 



where Ci — (cos(-\/37r/262)). It is beheved that for 
Kg > 1/4, the dSG model exhibits a second-order Ising 
transition by tuning the coupling constants. Indeed, as 
5 passes through a critical value 6c, the potential of the 
dSG model changes from a double-well type to a single- 
bottom one. We thus conclude that the transition at 
(5 = 5c is in the Ising universality class. The semiclassical 
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FIG. 4: Phase diagram of the spin-| tube. J_l(< J^) is fixed 
in the weak-rung-couping regime. The bosonization approach 
is less valid near both the lower and the upper critical fields. 
The upper critical field here is evaluated by the energy dis- 
persion of one magnon. 



analysis of Vs becomes more reliable with Vs being more 
relevant. It is hence expected that the value 5c mono- 
tonically decreases and approaches 6ci when Kgi^i K) 
increases— i.e., when H becomes close to the upper crit- 
ical field. Equation (3.11) and Fig. 12 in Ref. [12 also 
support this prediction. 

Since {1^3 j) oc (sin(\/27r6i)) « \/27r(6i) holds near 
the transition, the chirality would play the role of the 
Ising order parameter. Comparing our dSG model (fTU)) 
and results [Eqs. (29) and (58)] in Ref. [2l|, we can find 



i^h) {6 ^ 6cy^ 



(11) 



near 6 6c + 0. [The mean-field analysis of Vs leads 
to {1^3 j) ^ {6 — (5c)^/^.] The chirality correlation func- 
tion behaves as (k| jK^ 0) ~ C'2/ljl"'^^'* just at the critical 
point. Moreover, near the point, it follows the results in 
Ref. M- 

Sakai et al. have recently predicted that when 6 is 
added to the spin tube in zero field, the spin gap rapidly 
disappears Combining this prediction and ours, we 
construct the phase diagram on the 6-H plane as in 
Fig.H 



V. CONCLUSIONS AND DISCUSSIONS 

We have studied the three- leg AF spin-i tube ([T]) with 
a weak rung coupling under a magnetic field making use 
of Abelian bosonization techniques. In Sec. IIIH we have 
definitely elucidated without the help of any artificial 
approximation that the vector chiral long-range order 
occurs in the magnetic-field-driven one-component TLL 
phase. Due to the chiral order, the discrete rung-parity 
symmetry is spontaneously broken, but the U{1) spin- 
rotational symmetry remains unbroken. We have also 
been predicted that the chiral order survives even in the 
intermediate plateau state with M = 1/6. It is remark- 
able that one can obtain both gapless and gapped chiral 
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states by changing the strength of the magnetic field. In 
Sec. IIV[ we have investigated effects of the rung defor- 
mation in Eq. ([8|). It has been shown that when a posi- 
tive S is introduced with fixing J±d the chiral order and 
the TLL always continue, whereas when S is decreased 
the chiral order vanishes at a certain Ising critical point 
5 — 6c < and then a standard TLL phase appears (see 
Fig. [?]) . This implies that vector chirality correlations are 
favored in two-leg spin ladders^^ as well as in three-leg 
tubes. 

Quite recently, based on the spin- wave picture, we have 
predicted that the same chiral order still survives a cer- 
tain regime in the vicinity of saturation of the spin tube 
with 6 = 0^ Therefore, it is inferred that the tube 
possesses the chiral order in a quite wide range of the 
magnetic field. One immediately finds that the classical 
three-leg spin tube always has a vector chiral order from 
the zero field to the upper critical one. We thus may 
say that the classical nature is strong in the weak-rung- 
coupling area of the quantum spin tube Actually, 
as well known, the AF spin-i chain, which is the start- 
ing point of our analysis, has a large instability toward a 
classical Neel ordering. 

The discussion in this paper tells us that even without 
magnetic fields, a vector chiral order could also occur 
if a spin tube consists of three spin chains with a large 
K{> 1/2): for instance, the condition K > 1/2 is sat- 
isfied in a chain with easy-plane XX Z anisotropy. It is 



predicted that chiral orders exist in a certain parame- 
ter regime of easy-plane zigzag spin chains;^ Our results 
also suggests that a field-induced vector chiral order can 
appear in other frustrated odd-leg spin tubes— i.e., five-, 
seven-, nine-leg tubes, etc., tubes. 

Finally, we briefiy mention the possibility of experi- 
mentally detecting the vector chirality. The chiral order 
predicted in this paper must be destroyed by an effect 
of thermal fluctuation. However, chirality correlations 
are expected to still be strong if the temperature is suf- 
ficiently low. Moreover, weak three-dimensional inter- 
actions among tubes can stabilize the long-range chiral 
order: the gapless chiral order in the tube is expected to 
change into a conventional umbrella spin structure due 
to the interactions. Four-point spin correlation functions 
include features of the above chiral order or strong chiral 
correlation. Such functions, in principle, are measured in 
polarized neutron scattering, electromagnetic-wave reso- 
nance, etc. 
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